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1. INTRODUCTION 



■ Using the formalism of rigorous statistical mechanics, we study the phenom- 
^ I ena of phase separation and freezing-point depression upon freezing of solu- 
tions. Specifically, we devise an Ising-based model of a solvent-solute system 

I and show that, in the ensemble with a fixed amount of solute, a macroscopic 

■ phase separation occurs in an interval of values of the chemical potential of the 
I solvent. The boundaries of the phase separation domain in the phase diagram 

■ are characterized and shown to asymptotically agree with the formulas used in 

heuristic analyses of freezing point depression. The limit of infinitesimal con- 

• centrations is described in a subsequent paper 
Ph' 

> 

B 

> ! 1.1. Motivation 

• i— I , 

^ ■ The statistical mechanics of pure systems — most prominently the topic of 

^ ■ phase transitions and their associated surface phenomena — has been a sub- 

ject of fairly intensive research in recent years. Several physical principles for 
pure systems (the Gibbs phase rule, Wulff construction, etc.) have been put 
on a mathematically rigorous footing and, if necessary, supplemented with ap- 
propriate conditions ensuring their validity. The corresponding phenomena in 
systems with several mixed components, particularly solutions, have long been 
well-understood on the level of theoretical physics. However, they have not re- 
ceived much mathematically rigorous attention and in particular have not been 
derived rigorously starting from a local interaction. A natural task is to use the 
ideas from statistical mechanics of pure systems to develop a higher level of 
control for phase transitions in solutions. This is especially desirable in light 
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of the important role that basic physics of these systems plays in sciences, both 
general (chemistry, biology, oceanography) and applied (metallurgy, etc.). See 
e.g. [11,24,27] for more discussion. 

Among the perhaps most interesting aspects of phase transitions in mixed 
systems is a dramatic phase separation in solutions upon freezing (or boiling). 
A well-known example from "real world" is the formation of brine pockets in 
frozen sea water. Here, two important physical phenomena are observed: 

(1) Migration of nearly all the salt into whatever portion of ice/water mixture 
remains liquid. 

(2) Clear evidence of facetting at the water-ice boundaries. 
Quantitative analysis also reveals the following fact: 

(3) Salted water freezes at temperatures lower than the freezing point of pure 
water. This is the phenomenon of freezing point depression. 

Phenomenon (1) is what "drives" the physics of sea ice and is thus largely 
responsible for the variety of physical effects that have been observed, see 
e.g. [17, 18]. Notwithstanding, (1-3) are not special to the salt-water system; 
they are shared by a large class of the so called non-volatile solutions. A dis- 
cussion concerning the general aspects of freezing/boiling of solutions — often 
referred to as colligative properties — can be found in [24, 27]. 

Of course, on a heuristic level, the above phenomena are far from mys- 
terious. Indeed, (1) follows from the observation that, macroscopically, the 
liquid phase provides a more hospitable environment for salt than the solid 
phase. Then (3) results by noting that the migration of salt increases the en- 
tropic cost of freezing so the energy-entropy balance forces the transition point 
to a lower temperature. Finally, concerning observation (2) we note that, due to 
the crystalline nature of ice, the ice-water surface tension will be anisotropic. 
Therefore, to describe the shape of brine pockets, a Wulff construction has to 
be involved with the caveat that here the crystalline phase is on the outside. In 
summary, what is underlying these phenomena is a phase separation accom- 
panied by the emergence of a crystal shape. In the context of pure systems, 
such topics have been well understood at the level of theoretical physics for 
quite some time [12, 16, 32, 33] and, recently (as measured on the above time 
scale), also at the level of rigorous theorems in two [2, 4, 14, 22, 28, 29] and 
higher [6,9, 10] dimensions. 

The purpose of this and a subsequent paper is to study the qualitative 
nature of phenomena (1-3) using the formalism of equilibrium statistical me- 
chanics. Unfortunately, a microscopically realistic model of salted water/ice 
system is far beyond reach of rigorous methods. (In fact, even in pure water. 
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the phenomenon of freezing is so complex that crystaUzation in reahstic models 
has only recently — and only marginally — been exhibited in computer simula- 
tions [26].) Thus we will resort to a simplified version in which salt and both 
phases of water are represented by discrete random variables residing at sites 
of a regular lattice. For these models we show that phase separation dominates 
a non-trivial region of chemical potentials in the phase diagram — a situation 
quite unlike the pure system where phase separation can occur only at a single 
value (namely, the transition value) of the chemical potential. The boundary 
lines of the phase-separation region can be explicitly characterized and shown 
to agree with the approximate solutions of the corresponding problem in the 
physical-chemistry literature. 

The above constitutes the subject of the present paper. In a subsequent 
paper [1] we will demonstrate that, for infinitesimal salt concentrations scaling 
appropriately with the size of the system, phase separation may still occur dra- 
matically in the sense that a non-trivial fraction of the system suddenly melts 
(freezes) to form a pocket (crystal). In these circumstances the amount of salt 
needed is proportional to the boundary of the system which shows that the on- 
set of freezing-point depression is actually a surface phenomenon. On a qual- 
itative level, most of the aforementioned conclusions should apply to general 
non- volatile solutions under the conditions when the solvent freezes (or boils). 
Notwithstanding, throughout this and the subsequent paper we will adopt the 
language of salted water and refer to the solid phase of the solvent as ice, to 
the liquid phase as liquid-water, and to the solute as salt. 

1.2. General Hamlltonlan 

Our model will be defined on the rf-dimensional hypercubic lattice Z'*. We will 
take the (formal) nearest-neighbor Hamiltonian of the following form: 

pj^f = - ^{ajjy + a^L^Ly) + k ^ - jJ-sSx - tJ-J^x- (1-1) 

{x,y) XXX 

Here (5 is the inverse temperature (henceforth incorporated into the Hamito- 

nian), x and y are sites in 1/ and {x,y) denotes a neighboring pair of sites. 
The quantities 4, and are the ice (water), liquid (water) and salt vari- 
ables, which will take values in {0, 1} with the additional constraint 

4 + L, = l (1.2) 

valid at each site x. We will say that Ix — i indicates the presence of ice at x 
and, similarly, the presence of liquid at x. Since a single water molecule 
cannot physically be in an ice state, it is natural to interpret the phrase 4 = 1 
as referring to the collective behavior of many particles in the vicinity of x 
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which are enacting an ice-Uke state, though we do not formally incorporate 
such a viewpoint into our model. 

The various terms in (1.1) are essentially self-explanatory: An interac- 
tion between neighboring ice points, similarly for neighboring liquid points 
(we may assume these to be attractive), an energy penalty k for a simultane- 
ous presence of salt and ice at one point, and, finally, fugacity terms for salt 
and liquid. For simplicity (and tractability), there is no direct salt-salt interac- 
tion, except for the exclusion rule of at most one salt "particle" at each site. 
Additional terms which could have been included are superfluous due to the 
constraint (1.2). We will assume throughout that k > 0, so that the salt- ice 
interaction expresses the negative affinity of salt to the ice state of water. This 
term is entirely — and not subtly — responsible for the general phenomenon of 
freezing point depression. We remark that by suitably renaming the variables, 
the Hamiltonian in (1.1) would just as well describe a system with boiling point 
elevation. 

As we said, the variables 4 and indicate the presence of ice and liq- 
uid water at site x, respectively. The assumption + — 1 guarantees that 
something has to be present at x (the concentration of water in water is unity); 
what is perhaps unrealistic is the restriction of Ij. and to only the extreme 
values, namely Ix^L^ G {0, 1}. Suffice it to say that the authors are confident 
(e.g., on the basis of [3]) that virtually all the results in this note can be ex- 
tended to the cases of continuous variables. However, we will not make any 
such mathematical claims; much of this paper will rely heavily on preexisting 
technology which, strictly speaking, has only been made to work for the dis- 
crete case. A similar discussion applies, of course, to the salt variables. But 
here our restriction to e {0, 1} is mostly to ease the exposition; virtually 
all of our results directly extend to the cases when takes arbitrary (positive) 
real values according to some a priori distribution. 

1.3. Reduction to Ising variables 

It is not difficult to see that the "ice-liquid sector" of the general Hamiltonian 
(1.1) reduces to a ferromagnetic Ising spin system. On a formal level, this is 
achieved by passing to the Ising variables — — 4, which in light of the 
constraint (1.2) gives 

\ -\- Gx 1 — Cj; 

i^x = — ^ — and Ix = — - — . (1.3) 
By substituting these into (1.1), we arrive at the interaction Hamiltonian: 

I3M' = - J ^(TxCTy - h^ax + K^Sx^—^^ -^UsSx, (1.4) 

{x,y) XX X 
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where the new parameters J and h are given by 

J = — - — and h = -(ckl - aj) + —. (1.5) 

We remark that the third sum in (1.4) is still written in terms of "ice" indicators 
so that will have a well defined meaning even if k = oo, which corresponds 
to prohibiting salt entirely at ice-occupied sites. (Notwithstanding, the bulk of 
this paper is restricted to finite k.) Using an appropriate restriction to finite 
volumes, the above Hamitonian allows us to define the corresponding Gibbs 
measures. We postpone any relevant technicalities to Section 2.1. 

The Hamiltonian as written foretells the possibility of fluctuations in the 
salt concentration. However, this is not the situation which is of physical in- 
terest. Indeed, in an open system it is clear that the salt concentration will, 
eventually, adjust itself until the system exhibits a pure phase. On the level of 
the description provided by (1.4) it is noted that, as grand canonical variables, 
the salt particles can be explicitly integrated, the result being the Ising model 
at coupling constant J and external field h^^, where 

1 1 + e'*' 
Kn^h-r-\og ^^^^_^ . (1.6) 

In this context, phase coexistence is confined to the region /ieff = 0, i.e., a 
simple curve in the (/Xs, /i)-plane. Unfortunately, as is well known [5, 19, 20, 
23, 30], not much insight on the subject of phase separation is to be gained 
by studying the Ising magnet in an external field. Indeed, under (for example) 
minus boundary conditions, once h exceeds a particular value, a droplet will 
form which all but subsumes the allowed volume. The transitional value of h 
scales inversely with the linear size of the system; the exact constants and 
the subsequent behavior of the droplet depend on the details of the boundary 
conditions. 

The described "failure" of the grand canonical description indicates that 
the correct ensemble in this case is the one with a fixed amount of salt per unit 
volume. (The technical definition uses conditioning from the grand canonical 
measure; see Section 2.1.) This ensemble is physically more relevant because, 
at the moment of freezing, the salt typically does not have enough "mobility" 
to be gradually released from the system. It is noted that, once the total amount 
of salt is fixed, the chemical potential ji^ drops out of the problem — the relevant 
parameter is now the salt concentration. As will be seen in Section 2, in our 
Ising-based model of the solvent-solute system, fixing the salt concentration 
generically leads to sharp phase separation in the Ising configuration. More- 
over, this happens for an interval of values of the magnetic field h. Indeed, 
the interplay between the salt concentration and the actual external field will 
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demand a particular value of the magnetization, even under conditions which 
will force a droplet (or ice crystal, depending on the boundary condition) into 
the system. 

Remark 1.1. We finish by noting that, while the parameter h is for- 
mally unrelated to temperature, it does to a limited extent play the role of tem- 
perature in that it reflects the a priori amount of preference of the system for 
water v* ice. Thus the natural phase diagram to study is in the (c, /i) -plane. 

1 .4. Heuristic derivations and outline 

The reasoning which led to formula (1.6) allows for an immediate heuristic ex- 
planation of our principal results. The key simplification — which again boils 
down to the absence of salt-salt interaction — is that for any Ising configura- 
tion, the amalgamated contribution of salt, i.e., the Gibbs weight summed over 
salt configurations, depends only on the overall magnetization and not on the 
details of how the magnetization gets distributed about the system. In systems 
of linear scale L, let 3l(^) denote the canonical partition function for the 
Ising magnet with constrained overall magnetization M. The total partition 
function Zl{c, h) at fixed salt concentration c can then be written as 

Z^(c,/i) = ^3L(M)e'^^iyaM,c), (1.7) 

M 

where Wl{M, c) denotes the sum of the salt part of the Boltzmann weight — 
which only depends on the Ising spins via the total magnetization M — over all 
salt configurations with concentration c. 

As usual, the physical values of the magnetization are those bringing the 
dominant contribution to the sum in (1.7). Let us recapitulate the standard 
arguments by first considering the case c = (which implies Wl = 1), i.e., 
the usual Ising system at external field h. Here we recall that '^^{mL'^) can 
approximately be written as 

3L(mL'^) ~e-^'[-^''("^)+^l, (1.8) 

where C is a suitably chosen constant and ^j{m) is a (normalized) canonical 
free energy. The principal fact about ,'^j{m) is that it vanishes for m in the 
interval [— m^,m^], where = m^,{J) denotes the spontaneous magnetiza- 
tion of the Ising model at coupling J, while it is strictly positive and strictly 
convex for m with \m\ > m*. The presence of the "flat piece" on the graph 
of ^j{m) is directly responsible for the existence of the phase transition in the 
Ising model: For h > the dominant contribution to the grand canonical parti- 
tion function comes from M > m^L'^ while for /i < the dominant values of 



Colligative properties of solutions 



7 



the overall magnetization are M < —m^L'^. Thus, once — m^{J) > — 
which happens for J > Jc{d) with Jc{d) e (0, oo) whenever d > 2 — a phase 
transition occurs at h = 0. 

The presence of salt variables drastically changes the entire picture. In- 
deed, as we will see in Theorem 2.1, the salt partition function Wl{M,c) 
will exhibit a nontrivial exponential behavior which is characterized by 
a strictly convex free energy. The resulting exponential growth rate of 
'5L{M)e^^WL{M, c) for M ^ mL'^ is thus no longer a function with a flat 
piece — instead, for each h there is a unique value of ni that optimizes the 
corresponding free energy. Notwithstanding (again, due to the absence of 
salt-salt interactions) once that m has been selected, the spin configurations 
are the typical Ising configurations with overall magnetizations M ^ mL'^. 
In particular, whenever Zl{c^ h) is dominated by values of M mL'^ for 
an m G (— 'n^*, m^,), a macroscopic droplet develops in the system. Thus, due 
to the one-to-one correspondence between h and the optimal value of m, phase 
separation occurs for an interval of values of h at any positive concentration; 
see Fig. 1. 

We finish with an outline of the remainder of this paper and some discus- 
sion of the companion paper [1]. In Section 2 we define precisely the model 
of interest and state our main results concerning the asymptotic behavior of the 
corresponding measure on spin and salt configurations with fixed concentra- 
tion of salt. Along with the results comes a description of the phase diagram 
and a discussion of freezing-point depression, phase separation, etc., see Sec- 
tion 2.3. Our main results are proved in Section 3. In [1] we investigate the 
asymptotic of infinitesimal salt concentrations. Interestingly, we find that, in 
order to induce phase separation, the concentration has to scale at least as the 
inverse linear size of the system. 

2. RIGOROUS RESULTS 
2.1. The model 

With the (formal) Hamiltonian (1.4) in mind, we can now start on developing 
the mathematical layout of the problem. To define the model, we will need 
to restrict attention to finite subsets of the lattice. We will mostly focus on 
rectangular boxes C Z"' of L x L x ■ ■ • x L sites centered at the origin. 
Our convention for the boundary, (9A, of the set A C Z'' will be the collection 
of sites outside A with a neighbor inside A. For each x e A, we have the 
water and salt variables, ax G {— 1,+1} and Sx E {0,1}. On the boundary, we 
will consider fixed configurations ag^', most of the time we will be discussing 
the cases crgA = -1-1 or ct^a = —1, referred to as plus and minus boundary 
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conditions. Since there is no salt-salt interaction, we may as well set S^: — 

for all a; G 

We will start by defining the interaction Hamiltonian. Let A C Z'^ be a 
finite set. For a spin configuration aQA and the pair {a a, Sa) of spin and salt 
configurations, we let 

/3J^((7a, SA\agA) 



Here, as before, (,x, y) denotes a nearest-neighbor pair on Z'^ and the parame- 
ters J, h and k are as discussed above. (In light of the discussion from Sec- 
tion 1.3 the last term in (1.4) has been omitted.) The probability distribution of 
the pair (a a, Sa) takes the usual Gibbs-Boltzmann form: 

g-/3,i{*5\{o-A,SAkaA) 
^A(o'aA) 

where the normalization constant, Za{o-qa), is the partition function. The dis- 
tributions in Al with the plus and minus boundary conditions will be denoted 
by and , respectively. 

For reasons discussed before we will be interested in the problems with a 
fixed salt concentration c e [0, 1]. In finite volume, we take this to mean that 
the total amount of salt, 

Nl = A^l(S) = (2.3) 

is fixed. To simplify future discussions, we will adopt the convention that "con- 
centration c" means that < c\Al\ < Nl + 1, i.e., = [cL'^\. We may 
then define the finite volume Gibbs probability measure with salt concentra- 
tion c and plus (or minus) boundary conditions denoted by P^''^'^ (or P^''^''*). 
In light of (2.2), these are given by the formulas 

PL^'\-)-Pi{-\NL^[cL'\). (2.4) 

Both measures P^''^''* depend on the parameters J and k in the Hamiltonian. 
However, we will always regard these as fixed and suppress them from the 
notation whenever possible. 

2.2. Main theorems 

In order to describe our first set of results, we will need to bring to bear a few 
standard facts about the Ising model. For each spin configuration a = [gx) G 



-J 



{x,y) x&\ x&\ 



1 - (Tx 



(2.1) 
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{—1, 1}^^ let us define the overall magnetization in by the formula 

Ml = MM = (2.5) 

xeA^ 

Let m(/i, J) denote the magnetization of the Ising model with coupling con- 
stant J and external field /i > . As is well known, cf the proof of Theorem 3.1, 
/i I— > m(/i, J) continuously (and strictly) increases from the value of the sponta- 
neous magnetization = m(0, J) to one as h sweeps through [0, oo). In par- 
ticular, for each m G [m(0, J), 1), there exists a unique f) = I) (m, J) e [0, oo) 
such that m([) , J) = m. 

Next we will use the above quantities to define the func- 
tion ^j: (—1,1) — > [0,oo), which represents the canonical free energy 
of the Ising model in (1.8). As it turns out — see Theorem 3.1 in Section 3 — we 
simply have 

^j{m) = J dm'i){m',J)^{m^<m'<\m\}, me (-1,1). (2.6) 

As already mentioned, if J > Jc, where Jc = Jc{d) is the critical coupling con- 
stant of the Ising model, then > and thus (m) = for m G [— m*, m*] . 
(Since Jc{d) < oo only for d > 2, the resulting "flat piece" on the graph 
of m I— > ^j{m) appears only in dimensions d > 2.) From the perspective of 
the large-deviation theory, cf [13, 21], m i-^ ^j{m) is the large-deviation rate 
function for the magnetization in the (unconstrained) Ising model; see again 
Theorem 3.1. 

Let S^{p) = plogp + (1 — p) log(l — p) denote the entropy function of 
the Bernoulli distribution with parameter p. (We will set S^{p) — oo whenever 
p ^ [0, 1].) For each m e (-1, 1), each c e [0, 1] and each 9 e [0, 1], let 

^ ' ' ^ 2 Vl + m/ 2 V 1-m / 

As we will show in Section 3, this quantity represents the entropy of configu- 
rations with fixed salt concentration c, fixed overall magnetization m and fixed 
fraction 9 of the salt residing "on the plus spins" (and fraction 1 — 6 "on the 
minus spins"). 

Having defined all relevant quantities, we are ready to state our results. 
We begin with a large-deviation principle for the magnetization in the mea- 
sures Pl'"'^: 
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Theorem 2.1 . Let J > andn > be fixed. For each c e (0, 1), each 
/i e R and each m e (—1, 1), we have 

lim lim log Pf \ Ml - mL'^\ < eL^) 

^-Gn,c{m)+ inf Gh,c{m'). (2.8) 

m'G(— 1,1) 

Here m h- > Gh,c{p^) is given by 

Gh,c{m)- inf ^fe,e(m,^), (2.9) 
ee[o,i] 

e) = -/im - K^c - S(m, ^; c) + ^j{m). (2.10) 

The function m ^ Gh,c{fn) is finite and strictly convex on (—1,1) with 
lim^^-i-i G'y^ J^m) = ±oo. Furthermore, the unique minimizer m = m{h, c) of 
m ^ G/i_c(m) is continuous in both c and h and strictly increasing in h. 

On the basis of the above large-deviation result, we can now character- 
ize the typical configurations of the measures p^^'^'^. Consider the Ising model 
with coupling constant J and zero external field and let P^'"^ be the correspond- 
ing Gibbs measure in volume and ±-boundary condition. Our main result 
in this section is then as follows: 

Theorem 2.2. Let J > and k. > be fixed. Let c G (0,1) and h EM., 
and define two sequences of probability measures on [—1, 1] by the formula 

pt{[-lM) = PL"'\ML<mL''), me [-1,1]. (2.11) 

The measures allow us to write the spin marginal of the measure p^'*^'^ as 
a convex combination of the Ising measures with fixed magnetization; i.e., for 
any set A of configurations {(Tx)xeAL' have 

Pt'^^^A X {0, 1}^-) = j p±(dm) ¥1''{A\Ml = [mL^l). (2.12) 

Moreover, ifm = ■m{h, c) denotes the unique minimizer of the function m ^ 
Gh,c{'>TT')ffom (2.9), then the following properties are true: 

(1) Given the spin configuration on a finite set A C Z^, the salt variables 
on A are asymptotically independent. Explicitly, for each finite set A C 
Z'^ and any two configurations S\ e {0, 1}^ and a\ e {—1, 1}^, 

lim P±'^'^(SA = SA|fTA = frA) 

= + (1 - g^J5o(S.)}, (2.13) 
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where the numbers q± e [0,1] are uniquely determined by the equations 

(1+ Q- K . 1+m 1-m 

= e and q+ + g- = c. (2.14) 

1 — g+ 1 — g_ 2 2 



(2) 77je measure converges weakly to a point mass atm — m{h, c), 

lim = (2.15) 

/n particular, the Ising-spin marginal of the measure p^''^'^ is asymptot- 
ically supported on the usual Ising spin configurations with the overall 
magnetization Ml = (m + o{l))L'^, where m minimizes m i— > Gh^c.{'^)- 

The fact that conditioning p^ *^'^ on a fixed value of magnetization 
produces the Ising measure under same conditioning — which is the content 
of (2.12) — is directly related to the absence of salt-salt interaction. The princi- 
pal conclusions of the previous theorem are thus parts (1) and (2), which state 
that the presence of a particular amount of sdAi forces the Ising sector to choose 
a particular value of magnetization density. The underlying variational prin- 
ciple provides insight into the physical mechanism of phase separation upon 
freezing of solutions. (We refer the reader back to Section 1.4 for the physical 
basis of these considerations.) 

We will proceed by discussing the consequences of these results for the 
phase diagram of the model and, in particular, the phenomenon of freezing 
point depression. Theorems 2.1 and 2.2 are proved in Section 3.2. 



2.3. Phase diagram 

The representation (2.12) along with the asymptotic (2.15) allow us to char- 
acterize the distribution p^ '^ '* in terms of the canonical ensemble of the Ising 
ferromagnet. Indeed, these formulas imply that the distribution of Ising spins 
induced by P^''^'^ is very much like that in the measure P^'"^ conditioned on 
the event that the overall magnetization Ml is near the value m{h,c)L'^. As a 
consequence, the asymptotic statements (e.g., the Wulff construction) that have 
been (or will be) established for the spin configurations in the Ising model with 
fixed magnetization will automatically hold for the spin-marginal of the P^'"^''* 
as well. 

A particular question of interest in this paper is phase separation. Recall 
that = J) denotes the spontaneous magnetization of the Ising model 
at coupling J. Then we can anticipate the following conclusions about typical 
configurations in measure P^''^''*: 



12 



Alexander, Biskup and Chayes 




Fig. 1. The phase diagram of the ice- water system with k ^ 1. The horizontal axis marks 
the concentration of the salt in the system, the vertical line represents the external field acting 
on the Ising spins — see formula (1.5). For positive concentrations c > 0, the system stays in 
the liquid-water phase throughout a non-trivial range of negative values of h — a manifestation 
of the freezing-point depression. For {h, c) in the shaded region, a non-trivial fraction of the 
system is frozen into ice. Once {h, c) is on the left of the shaded region, the entire system is in 
the ice state. For moderate values of k, the type of convexity of the transition lines may change 
from concave to convex near {h, c) = (0, 0); see the companion paper [1]. 

(1) If m(/i, c) > m*, then the entire system (with plus boundary condition) 
will look like the plus state of the Ising model whose external field is 
adjusted so that the overall magnetization on the scale is roughly 

m{h, c)L'^. 

(2) If m{h, c) < — m*, then the system (with minus boundary condition) 
will look like the Ising minus state with similarly adjusted external field. 

(3) If m{h,c) e (— m*,m*), then, necessarily, the system exhibits phase 
separation in the sense that typical configurations feature a large droplet 
of one phase inside the other. The volume fraction taken by the droplet 
is such that the overall magnetization is near m(/i, c)L'^. The outer phase 
of the droplet agrees with the boundary condition. 

The cases (1-2) with opposite boundary conditions — that is, the minus bound- 
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ary conditions in (1) and the plus boundary conditions in (2) — are still as stated; 
the difference is that now there has to be a large contour near the boundary flip- 
ping to the "correct" boundary condition. 

Remark 2.3. We have no doubt that the aforementioned conclusions 

(1-3) hold for all > 2 and all J > (with a proper definition of the droplet 
in part (3), of course). However, the depth of conclusion (3) depends on the 
level of understanding Wulff construction, which is at present rather different in 
dimensions d — 2 and d > 3. Specifically, while 'md — 2 the results of [14,22] 
allow us to claim that for all J > and all magnetizations m e (— m^, nrii), 
the system will exhibit a unique large contour with appropriate properties, in 
c? > 3 this statement is known to hold [6, 10] only in "L^-sense" and only 
for m e (— m*, m^) which are near the endpoints. (Moreover, not all values 
of J > Jc are, in principle, permitted; cf [7] for a recent improvement of these 
restrictions.) We refer to [8] for an overview of the situation. 

Notwithstanding the technical difficulties of Wulff construction, the above 
allows us to characterize the phase diagram of the model at hand. As indi- 
cated in Fig. 1, the /i < and c > quadrant splits into three distinct parts: 
The liquid-water region, the ice region and the phase separation region, which 
correspond to the situations in (1-3), respectively. The boundary lines of the 
phase-separation region are found by setting 

m(/i, c) = ±m*, (2.16) 

which in light of strict monotonicity oi h ^ m{h, c) allows us to calculate h 
as a function of c. The solutions of (2.16) can be obtained on the basis of the 
following observation: 

Proposition 2.4. Let m G [— m^,m^] and c G [0, 1] and define the 
quantities q± = g±(m, c, k) by (2.14). Let h be the solution to m{h, c) = m. 
Then we have: 

/i^ilogi^. (2.17) 
2 1 — 

In particular, there exist two continuous and decreasing functions 
h±: [0, oo) — > (— cx),0] with /i+(c) > /i-(c) for all c > 0, such that 
— < m{h, c) < is equivalent to h-{c) < h < h+{c) for all c > 0. 

Proposition 2.4 is proved at the very end of Section 3.2. Here is an infor- 
mal interpretation of this result: The quantities q± represent the mole fractions 
of salt in liquid-water and ice, respectively. In mathematical terms, g+ is the 
probability of having a salt particle on a given plus spin, and g_ is the corre- 
sponding quantity for minus spins, see (2.13). Formula (2.17) quantifies the 
shift of the chemical potential of the solvent (which is given by 2h in this case) 
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due to the presence of the solute. This is a manifestation of freezing point 
depression, see also Remark 1.1. In the asymptotic when c <^ 1 we have 

2hf^q--q+. (2.18) 

This relation, derived in standard chemistry and physics books under the aus- 
picies of the "usual approximations," is an essential ingredient in the classical 
analyses of coUigative properties of solutions [24,27]. Here the derivation is 
a direct consequence of a microscopic (albeit simplistic) model which further 
offers the possibility of systematically calculating higher-order corrections. 

3. PROOFS 

The proofs of our main results are, more or less, straightforward exercises in 
large-deviation analysis of Gibbs distributions. We first state and prove a cou- 
ple of technical lemmas; the actual proofs come in Section 3.2. 

3.1. Preliminaries 

The starting point of the proof of Theorem 2.1 (and, consequently. Theo- 
rem 2.2) is the following large-deviation principle for the Ising model at zero 
external field: 

Tlieorem 3.1 . Consider the Ising model with coupling constant J G 
[0, oo) and zero external field. LetW'^''^ be the corresponding (grand canonical) 
Gibbs measure in volume A^, and -^-boundary conditions. Then for all m e 

[-1,1]. 

lim hm logP^'-^dMi - mL'^\ < eL'^) = -=Fj(m), (3.1) 

where Ml is as in (2.5) and is as defined in (2.6). 

Proof. The claim is considered standard, see e.g. [31, Section II.l], and 
follows by a straightforward application of the thermodynamic relations be- 
tween the free energy, magnetization and external field. For completeness (and 
reader's convenience) we will provide a proof. 

Consider the function ^^(/i) = logE^"^(e''^^), where E^"^ is the ex- 
pectation with respect to P^'"^, and let 4>{h) = limi^^ 4>L{h). The limit ex- 
ists by subadditivity arguments and is independent of the boundary condition. 
The function h ^ (j){h) is convex on R, real analytic (by the Lee- Yang theo- 
rem [25]) on R \ {0}, and hence it is strictly convex on R. By the h ^ —h 
symmetry there is a cusp at /i = whenever = (p'{0~^) > 0. It follows that 
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for each m e [m*, 1) there is a unique f) = f)(m, J) > such that — m, 
with [)(m, J) increasing continuously from to cxd as m increases from 
to 1. The plus-minus symmetry shows that a similar statement holds for the 
magnetizations in (—1, — m^]. 

Let (j)* denote the Legendre transform of 0, i.e., (j)*{m) = sup i^^^[mh — 
(f){h)]. By the above properties of /i i— > (f){h) we infer that 0*(m) = mi) — 0(f)) 
when m e (—1, — m*) U (m*, 1) and () = f)(m, J), while (t)*{m) = —0(0) = 
for m G [— m^,m^]. Applying the Gartner- Ellis theorem (see [21, Theo- 
rem V.6] or [13, Theorem 2.3.6]), we then have (3.1) with ,^j{m) = (f)*{m) 
for all m G [—1, — m^) U (m*, 1] — which is the set of so called exposed points 
of 0*. Since 0*(±m*) ~ and the derivative of m i— > (f)*{m) is i){m, J), 
this is given by the integral in (2.6). To prove (3.1) when m G [— m*, m^], 
we must note that the left-hand side of (3.1) is nonpositive and concave in m. 
(This follows by partitioning Al into two parts with their own private magneti- 
zations and disregarding the interaction through the boundary.) Since -^j(m) 
tends to zero as m tends to ±m^ we thus have that (3.1) for m G [ m^, m^] as 
well. I 

Remark 3.2. The "first" part of the Gartner-Ellis theorem [21, Theo- 
rem V.6] actually guarantees the following large-deviation principle: 

lim sup ^ log P^'-^ eC) <- inf 0*(m) (3.2) 

for any closed set C C R while 

liminf ^logP^'-'f^ G > - inf 0*(m) (3.3) 

for any open set O C R. (Here (f)*{m) = ^j{m) for m G [—1, 1] and (j)*{m) = 
oo otherwise.) The above proof follows by specializing to e-neighborhoods of 
a given m and letting e | 0. The m G [— my^, rrii] cases — i.e, the non-exposed 
points — have to be dealt with separately. 

The above is the core of our proof of Theorem 2.1. The next step will be 
to bring the quantities c and h into play. This, as we shall see, is easily done if 
we condition on the total magnetization. (The cost of this conditioning will be 
estimated by (3.1).) Indeed, as a result of the absence of salt-salt interaction, 
the conditional measure can be rather precisely characterized. Let us recall the 
definition of the quantity Nl from (2.3) which represents the total amount of 
salt in the system. For any spin configuration a = (a^) G {—1, 1}^^ and any 
salt configuration S = (S^) G (0, 1}^^, let us introduce the quantity 

Ql = Ql{<j, s)^J2 ^-^-^ ^^-"^^ 
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representing the total amount of salt "on the plus spins." Then we have: 

Lemma 3.3. For any fixed spin configuration a = (a^) G { — 1, 1}^^, 
all salt configurations (Sx) E {0, 1}^^ with the same Nl and Ql have the 
same probability in the conditional measure P^''^''^{-\a = a). Moreover, for 
any S = (S^) e {0, 1}"^^ with Nl = [cL'^\ and for any m e [-1, 1], 

Pj'''^(S occurs, Ml - [mL'^J) 

- ^{ML{a)=lmLdi}), P-Jj 

where the normalization constant is given by 

^L- 2^ l{^^(s')=LcL'ij}% ie '')■ P-oj 

Here E^''^ is the expectation with respect to P^''^. 

Proof. The fact that all salt configurations with given A^^, and Ql have 
the same probability in P^''^''^(- jcr = ct) is a consequence of the observation 
that the salt-dependent part of the Hamiltonian (2.1) depends only on Ql- The 
relations (3.5-3.6) follow by a straightforward rewrite of the overall Boltzmann 
weight. I 

The characterization of the conditional measure 
from Lemma 3.3 allows us to explicitly evaluate the configurational entropy 
carried by the salt. Specifically, given a spin configuration a — (ax) e 
{-1, 1}"^^ and numbers 0,ce (0, 1), let 

A'£'{a) = {{Sx) G {0, 1}^- : Nl = [cL"], Ql = L^cL^J }. (3.7) 

The salt entropy is then the rate of exponential growth of the size of .A^''^((t) 
which can be related to the quantity S(m, 9; c) from (2.7) as follows: 

Lemma 3.4. For each e' > and each rj > there exists a num- 
ber Lq < CO such that the following is true for any 9,c E (0, 1), any m e 
(— 1, 1) that obey \m\ < 1 — rj, 

* <1_, an, ?(i^<l-,, (3.8) 



1 + m 1 — m 

and any L > Lq: If a = (ax) G { — 1,1}^^ is a spin configuration with 
Ml{(t) = [mL'^J, then 



< e'. (3.9) 
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Proof. We want to distribute A^^ = L^^*^] ^^^t particles over posi- 
tions, such that exactly Ql = [9cL^\ of them land on ^{L'^ + M^) plus sites 
and Nl — Ql on ^{L^ — Ml) minus sites. This can be done in 

number of ways. Now all quantities scale proportionally to L'^ which, applying 
Stirling's formula, shows that the first term is within, say, e^^ ^'^^ multiples of 

once L > Lq, with Lq depending only on e'. A similar argument holds also for 
the second term with 6 replaced by 1 — and m by —m. Combining these ex- 
pressions we get that |w4^''^(cr) | is within e^^ multiples of exp{L'^S(m, Q\ c)} 
once L is sufficiently large. | 

For the proof of Theorem 2.2, we will also need an estimate on how many 
salt configurations in ^^{a) take given values in a finite subset A c A^. To 
that extent, for each a e {—1, 1}^^ and each Sa e {0, 1}^ we will define the 
quantity 

As a moment's thought reveals, R^/^];^{a, Sa) can be interpreted as the prob- 
ability that {Sa = Sa} occurs in (essentially) any homogeneous product 
measure on S = (S^) G {0, 1}"^^ conditioned to have Nl{S) = [cL'^\ 
and QLio', S) = l9cL^\. It is therefore not surprising that, for spin configu- 
rations (7 with given magnetization, -Ra ') will tend to a product measure 
on Sa e {0, 1}"^. A precise characterization of this limit is as follows: 

Lemma 3.5. For each e > 0, each K > 1 and each r) > there 
exists Lq < oo such that the following holds for all L > Lq, all A C Al 
with \ A\ < K, all m with \m\ < 1 — r] and all 9,c & [r),l — vlfar which 

29c 2{1 -9)c 

p+ = — and P- = — ; (3.13) 

1 + m 1 — m 

satisfy p± G [rj,! — r]]: If a = (a^;) G { — 1, 1}'^^ is a spin configuration such 
that Ml{cf) — [mL'^J and Sa G {0, 1}^ is a salt configuration in A, then 



< e. (3.14) 
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Proof. We will expand on the argument from Lemma 3.4. Indeed, 
from (3.10) we have an expression for the denominator in (3.12). As to the 
numerator, introducing the quantities 

M^ = Y,(^x, A^A = J]S^, Qa = J2s^^-^, (3.15) 

xeA xeA xeA 

and the shorthand 

r-£\ (r' - e 



D = Dr,r>,sA^, g, q') = l(Y,^ , (3.16) 



s 

the same reasoning as we used to prove (3.10) allows us to write the ob- 
ject R^^\{a^SfC) as Dry^s,s'{^,^',q,q')^ where the various parameters are as 
follows: The quantities 

r = and r = (3.17) 

2 2 

represent the total number of pluses and minuses in the system, respectively, 

s = Ql and s' = Nl-Ql (3.18) 
are the numbers of salt particles on pluses and minuses, and, finally, 

iJ^^Ml^, ^^J3^^ q^Q^ and q'^N^-Q^ (3.19) 

are the corresponding quantities for the volume A, respectively. 

Since (3.13) and the restrictions on |m| < 1 — rj and 9,c E [77, 1 — r^] imply 
that r, r', s, s', r — s and r' — s' all scale proportionally to L^, uniformly in a 
and Sa, while £ and f/ are bounded by |A| — which by our assumption is less 
than K — we are in a regime where it makes sense to seek an asymptotic form 
of quantity D. Using the bounds 

a'-e-^V" < i^±^ < a'e^y^ (3.20) 
a! 

which are valid for all integers a and b with |6| < a, we easily find that 

(1-;) {-) (l") +"(1). (3.21) 

Since s/r ^ p+ and s' /r' — > p_ as L — > 00, while £, q, £' and q' stay bounded, 
the desired claim follows by taking L sufficiently large. | 

The reader may have noticed that, in most of our previous arguments, 9 
and m were restricted to be away from the boundary values. To control the 
situation near the boundary values, we have to prove the following claim: 
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Lemma 3.6. For each e e (0, 1) and each L > 1, let be the event 

^L,e={|Mi|<(l-e)L^} 

n {e\{L'' + Ml) <Ql<{1- e)\{L'' + Ml)]. (3.22) 

Then for each c e (0, 1) and each /i e R there exists an e > such that 



1 



limsup ^ \ogP^'''^\£l^,) < 0. (3.23) 



L—>oo 



Proof. We will split the complement of £L,e into four events and prove 
the corresponding estimate for each of them. We begin with the event {Ml < 
— (1 — e)L'^}. The main tool will be stochastic domination by a product 
measure. Consider the usual partial order on spin configurations defined by 
putting a ^ a' whenever < a'^ for all x. Let 

A = inf min min Pt'^'^{a^ = S) (3.24) 

Se{-l,l}Ai 

be the conditional probability that cTj; = +1 occurs given a spin configuration a' 

in Kl \ {x} and a salt configuration S in Ax,, optimized over all a', S and 
also X ^ Kl and the system size. Since P^''^'^{o-x = 1\<7', S) reduces to (the 
exponential of) the local interaction between and its ultimate neighborhood, 
we have A > 0. 

Using standard arguments it now follows that the spin marginal of P^'"^''* 
stochastically dominates the product measure defined by Pa(cx = 1) = A 
for all x. In particular, we have 

P^'^'''{Ml < -(1 - e)L'^) < Pa (Ml < -(1 - e)L'^). (3.25) 

Let e < 2 A. Then A — (1 — A) — namely, the expectation of ax with respect 
to Pa — exceeds the negative of (1 — e) and so Cramer's theorem (see [21, 
Theorem L4] or [ 1 3, Theorem 2. 1 .24]) implies that the probability on the right- 
hand side decays to zero exponentially in L*^, i.e., 

limsup^ log Pa (Mx, < -(1 - e)L'=') < 0. (3.26) 

The opposite side of the interval of magnetizations, namely, the event {Ml > 
(1 — e)L'^}, is handled analogously (with A now focusing on ax — —1 instead 
of ax = 1). 

The remaining two events, marking when is either less than e or larger 
than (1 — e) times the total number of plus spins, are handled using a similar 
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argument combined with standard convexity estimates. Let us consider the 
event {Ql < eL"^}— which contains the event {Ql < e\{ML + L'^)}— and let 
us emphasize the dependence of the underlying probability distribution on k by 
writing P^''^'^ as P^. Let denote the expectation with respect to and note 
that E^(/) = Eo(/e'^^^)/Eo(e'^'3L), We begin by using the Chemoff bound to 
get 

^.{Ql < eL") < e«^^ E«(e-»«-) = — — , a > 0. (3.27) 

A routine application of Jensen's inequality gives us 

K(Ql < eL") < exp{a(eL'^ - E«_„(Q^))}. (3.28) 

It thus suffices to prove that there exists ?l k' < k such that infi>i j^d^K'iQL) 
is positive. (Indeed, we take e to be strictly less than this number and set a — 
K — k' X.O observe that the right-hand side decays exponentially in L'^.) To show 
this we write ^f^^Qi) as the sum of ^^'{(^x = 1, S^; = 1) over all x G Kl- 
Looking back at (3.24), we then have P«'((J^ = 1, = 1) > AP„/(S^ = 1), 
where A is now evaluated for k' , and so 

^n'{.QL) > A ^ f^,{S, = 1) = AE«,(7Vi) « AcL'^. (3.29) 

Thus, once Ac > e, the probability Pk((5l < ^L'^) decays exponentially in L'^. 

As to the complementary event, {Ql > ( 1 — e) \{M l + L^)}, we note that 
this is contained in {H^ < eL'^}, where Hl counts the number of plus spins 
with no salt on it. Since we still have E„(/) = Eo{fe-''^i^)/Eo{e-''^^), the 
proof boils down to the same argument as before. | 

3.2. Proofs of Theorems 2.1 and 2.2 

On the basis of the above observations, the proofs of our main theorems are 
easily concluded. However, instead of Theorem 2. 1 we will prove a slightly 
stronger result of which the large-deviation part of Theorem 2. 1 is an easy 
corollary. 

Theorem 3.7. Let J > and k, > be fixed. For each c, G (0, 1), 
each /?. G M and each m G (—1, 1), let Bl^^ = BL^^{m, c, 9) be the set of all 
(a, S) G {-1, 1}^^ X {0, Ij^^forwhich \ML-mL'^\ < tL'^ and \QL-ecL'^\ < 
eL'^ hold. Then 

lim lim \ = -%,c{m, 0) + inf %,c{m' , 9'), (3.30) 

ej.0 L— »oo iv" m'e(— 1,1) 

O'e[0,l] 
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where ^h,c("^) ^■s' in (2.10). 



Proof. Since the size of the set A^£'^{a) depends only on the overall mag- 



netization, let A^'' (m) denote this size for the configurations a with ML{a) — 
[mL^l . First we note that, by Lemma 3.3, 

P^'^'^Ql = iecL'\, M, = [mL'\) = El^^ (3.31) 

where 

Xi(m, Q) = A\\m) e'»LmL''J+«LecL<'J p±-^(Mi = [mL'^J). (3.32) 

Here Zj^ is the normalization constant from (3.6) which in the present formula- 
tion can also be interpreted as the sum of Ki{m, 9) over the relevant (discrete) 
values of m and 9. 

Let KL.c{'m,9) denote the sum of KL{m',9) over all m' and 9' for 
which m'L'^ and 9'cL'^ are integers and \m' — m\ < e and \9'c — 9c\ < e. 
(This is exactly the set of magnetizations and spin-salt overlaps contributing 
to the set Bl^^.) Applying (3.1) to extract the exponential behavior of the last 
probability in (3.32), and using (3.9) to do the same for the quantity ^^''^(m), 
we get 



<e + e', (3.33) 



where e' is as in (3.9). As a consequence of the above estimate we have 

,to,to!28%i!^ = (3.34) 

for any m G (—1, 1) and any 9 E (0, 1). 

Next we will attend to the denominator in (3.31). Pick 6 > and consider 
the set 

Ms ^ {{m,9): \m\ < 1 - 5, 5 < 9 < 1 - S}. (3.35) 

We will write Zl as a sum of two terms, Zl = Z^l^ + Z^j^\ with Z^j^^ oh- 
tained by summing K{m, 9) over the admissible (m, 9) G M.^ and Z\^ col- 
lecting the remaining terms. By Lemma 3.6 we know that zf^ /Zl decays 
exponentially in L'^ and so the decisive contribution to Zl comes from Z^l^ . 
Assuming that e <^ 5, let us cover by finite number of sets of the 
form [m^ — e, + e] x [6^^ — e,9'^ + e], where m'^ and 9'^ are such that m'^L'^ 
and 9'^cL'^ are integers. Then Z^^^ can be bounded as in 

maxXi,,(m^, 9',) < Zf < J] Xi,,(m^, 9',). (3.36) 
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Moreover, the right-hand side is bounded by the left-hand side times a poly- 
nomial in L. Taking logarithms, dividing by L^, taking the limit L oo, 
refining the cover and applying the continuity of (m, 9) ^hfii'm, allows 
us to conclude that 

lim-^^|-^ = - inf mi % Am, 9). (3.37) 

Combining these observations, (2.8) is proved. | 

Prooi of Theorem 2. 1. The conclusion (2.8) follows from (3.30) by 
similar arguments that prove (3.37). The only remaining thing to prove is 
the strict convexity of m i— > Gh,cim) and continuity and monotonicity of its 
minimizer. First we note that 9 ^i^cifn, 9) is strictly convex on the set 
of 9 where it is finite, which is a simple consequence of the strict convexity 
of p I— > y{p). Hence, for each m, there is a unique 9 — 9{m) which mini- 
mizes 9 I— > ^h^ci'lTT': 

Our next goal is to show that, for kc > 0, the solution 9 — 9{m) will 
satisfy the inequality 

> (3.38) 

(A heuristic reason for this is that 9 = corresponds to the situation when 
the salt is distributed independently of the underlying spins. This is the dom- 
inating strategy for k = 0; once k > it is clear that the fraction of salt 
on plus spins must increase.) A formal proof runs as follows: We first note 
that 9{m) solves for 9 from the equation 

d 

— S(m, 9] c) = -KC, (3.39) 

where 5(m, 9] c) is as in (2.7). But 9 i— > S(m, 9] c) is strictly concave and its 
derivative vanishes at ^ = \{1 -\-m). Therefore, for > the solution 9 — 
9{m) of (3.39) must obey (3.38). 

Let V be the set of (m, ^) G (-1, 1) x (0, 1) for which (3.38) holds and 
note that V is convex. A standard second-derivative calculation now shows 
that '^h^df^i 6') is strictly convex on V. (Here we actually differentiate the func- 
tion '^h,c{fn-, — =^j(?77.) — which is twice differentiable on the set where it is 
finite — and then use the known convexity of ^j(m). The strict convexity is 
violated on the line 9 = \{\ + m) where (m, 9) ^ ^h,c{'m, 9) has a flat piece 
for m e [— m*, m*].) Now, since 9{m) minimizes '^h,c{'m, 9) for a given m, the 
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strict convexity of ^/i_c(^; ^) on V implies that for any A e (0, 1), 

Gh,c{>^mi + (1 - A)m2) 

< $fv(Ami + (1 - A)m2, A^(mi) + (1 - A)^(m2)) 

< X%,c{mi, e{m,)) + (1 - A)^;,,e(m2, ^K)) (3.40) 
= AG/,,e(mi) + (1 - X)Gh,c{m2). 

Hence, m i-^ Gh,ci'm) is also strictly convex. The fact that G'{m) diverges 
as m — > ±1 is a consequence of the corresponding property of the func- 
tion m ^ ^j{m) and the fact that the rest of c is convex in m. 

As a consequence of strict convexity and the abovementioned "steepness" 
at the boundary of the interval (—1,1), the function m i— > Gh,c{'ni') has a unique 
minimizer for each h eM. and c > 0, as long as the quantities from (3.13) sat- 
isfy p± < 1. The minimizer is automatically continuous in h and is manifestly 
non-decreasing. Furthermore, the continuity of Gh,c in c allows us to conclude 
that 9{m) is also continuous in c. What is left of the claims is the strict mono- 
tonicity of m as a function of h. Writing Gh,c{Tn) as —hm + g{m) and noting 
that g is continuously differentiable on (—1, 1), the minimizing m satisfies 

g'(m) = h. (3.41) 

But g{m) is also strictly convex and so g'{m) is strictly increasing. It follows 
that m has to be strictly increasing with h. | 

Theorem 3.1 has the following simple consequence that is worth high- 
Ughting: 

Corollary 3.8. For given /i e M and c e (0, 1), let (m, 9) be the mini- 
mizer of^h^cii^-i d)- Then for all e > 0, 

lim Pl'^'^^HQl - 9cL'^\ > eL'^ or \Ml - mL'^\ > eL'^) = 0. (3.42) 

Proof. On the basis of (3.30) and the fact that '^h,c{m, 9) has a unique 
minimizer, a covering argument — same as used to prove (3.37) — implies that 
the probability on the left-hand side decays to zero exponentially in L'^. | 

Before we proceed to the proof of our second main theorem, let us make 
an observation concerning the values of p± at the minimizing m and 9: 

Lemma 3.9. Let /i G M and c G (0, 1) fixed and let (m, 9) be 
the minimizer of'^h,c{'^-i 9). Define the quantities q± = q±{m, c, n) by (2.14) 
andp± — p±{m, 9, c) by (3.13). Then 

q+ = p+ and = p^. (3.43) 

Moreover, q± are then related to h via (2.17) whenever m G [— m*, m*]. 
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Proof. First let us ascertain that q± are well defined from equations 
(2.14). We begin by noting that the set of possible values of q^) is the unit 
square [0, 1]^. As is easily shown, the first equation in (2.14) corresponds to an 
increasing curve in [0, 1]^ connecting the comers (0, 0) and (1,1). On the other 
hand, the second equation in (2. 14) is a straight line with negative slope which 
by the fact that c < 1 intersects both the top and the right side of the square. 
It follows that these curves intersect at a single point — the unique solution of 
(2.14). 

Next we will derive equations that p± have to satisfy. Let (m, 9) be the 
unique minimizer of '^h,c{^-, The partial derivative with respect to 9 yields 

c(j^'(p+) - S^'ip-)) = KC (3.44) 

and from the very definition of p± we have 

1 + m 1 — m ,^ .r-. 

P+ + -^P- = c. (3.45) 

Noting that y{p) = log we now see that p± satisfies the same equations 
as q± and so, by the above uniqueness argument, (3.43) must hold. 

To prove relation (2.17), let us also consider the derivative of '^h,c{'^i ^) 
with respect to m. For solutions in [— m*, m*] we can disregard the part 
of the function (because its vanishes along with its derivative throughout this 
interval), so we have 

/i=-^S(m,^;c). (3.46) 

A straightforward calculation then yields (2.17). | 
Now we are ready to prove our second main result: 

Proof of Thieorem 2.2. The crucial technical step for the present proof 
has already been established in Lemma 3.3. In order to plug into the latter 
result, let us note that the sum of e'**^^^'^''^) over all salt configurations S — 
(S-r) e {0, 1}^^ with Nl — [cL*^] is a number depending only on the total 
magnetization Ml — ML{a). Lemma 3.3 then implies 

p±'^''^(^ X {0, 1}^^ n {Ml = [mL''\}) 

^u;L{m)F^'\An{ML= ImL^}}) (3.47) 



where ujL{m) is a positive number depending on m, the parameters c, h, J 
and the bounc 
ply the distrit 
proves (2.12). 



and the boundary condition ± but not on the event A. Noting that is sim- 
ply the distribution of the random variables Ml/L'^ in measure this 
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In order to prove the assertion (2.13), we let a e {0, 1}^^, pick A c 
and fix S e {0, 1}^. Since Lemma 3.3 guarantees that, given {a — a}, all salt 
configurations with fixed Ql and concentration c have the same probability 

in P^''''^{-\a = a), we have 

P±'^''^(Sa = Sa, S e A^l%a)\a = a) = Sa), (3.48) 

where i?Ai is defined in (3.12). Pick 77 > and assume, as in Lemma 3.5, 
that c G [t],! — Tj], 6 E [r]A — rj] and Ml^o-) = \mL'^\ for some m with \m\ < 
1—1]. Then the aforementioned lemma tells us that 7?a'l(^' ') within e of 
the probability that Sa occurs in the product measure where the probability 
of S-r = 1 is p+ if ax — +1 and p_ if — —1. 

Let (m, 9) be the unique minimizer of '^h,c{fn, Taking expectation of 
(3.48) over a with (Ta fixed, using Corollary 3.8 to discard the events \Ml/L'^ — 
m\ > e or \Ql/ L''' — 6c\ > e and invoking the continuity of p± in m and 0, we 
find out that P^''^'^{S\ = S\\a\ = a\) indeed converges to 

Uip^JiiS.) + {l-p,MSx)}, (3.49) 

with p-t evaluated at the minimizing (m, 9). But for this choice Lemma 3.9 
guarantees thatp± = q±, which finally proves (2.13-2.14). | 

The last item to be proved is Proposition 2.4 establishing the basic features 
of the phase diagram of the model under consideration: 

Proof of Proposition 2.4. From Lemma 3.9 we already know that the 
set of points m{h, c) = m for m e [— m*, m*] is given by the equation (2.17). 
By the fact that m{h,c) is strictly increasing in h and that m{h,c) — > ±1 
as h ^ ±00 we thus know that (2. 17) defines a line in the (h, c)-plane. Special- 
izing to m = ±m^ gives us two curves parametrized by functions c ^^ h±{c) 
such that at (/i, c) satisfying h-{c) < h < h+{c) the system magnetiza- 
tion m{h, c) is strictly between — m^^ and m*, i.e., (/i, c) is in the phase sep- 
aration region. 

It remains to show that the above functions c ^ h±{c) are strictly mono- 
tone and negative for c > 0. We will invoke the expression (2. 17) which applies 
because on the above curves we have m{h, c) e [— m*, rrii]. Let us introduce 
new variables 

i?+ = and R_ = (3.50) 

1 - g+ 1 - g_ 

and, writing h in (2.17) in terms of R±, let us differentiate with respect to c. 
(We will denote the corresponding derivatives by superscript prime.) Since 
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(2.14) gives us that — e '^R+, we easily derive 

R'_ R' 1 — e"** 

^ TT^ " TT^ ^ "^+(i + i?+)(i + i?_)- ^^-^^^ 

Thus, /i' and i?^ have opposite signs; i.e., we want to prove that R'j^ > 0. 
But that is immediate: By the second equation in (2.14) we conclude that 
at least one of must be strictly positive, and by R- — e~'^R+ we find 
that both > 0. It follows that c i— > h±{c) are strictly decreasing, and 
since h±{0) = 0, they are also negative once c > 0. | 
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